Statistical fluctuations of ground-state energies 
and binding energies in nuclei 

A. Molinari and H. A. Weidenmiiller 
February 5, 2008 

Abstract 

The statistical fluctuations of the ground-state energy and of the bind- 
ing energy of nuclei are investigated using both perturbation theory and 
supersymmetry. The fluctuations are induced by the experimentally ob- 
served stochastic behavior of levels in the vicinity of neutron threshold. 
The results are compared with a recent analysis of binding-energy fluc- 
tuations by Bohigas and Leboeuf, and with theoretical work by Feshbach 
et al. 

Introduction. In this letter, we address the question: How does the inter- 
action with high-lying configurations affect the properties of nuclear ground 
states? In particular, what are the ensuing uncertainties, given the fact that 
not much is known about both, the interaction and the high-lying configura- 
tions? The question has arisen within two seemingly different contexts, (i) In a 
purely theoretical framework, and inspired by the analogy with the theory of the 
optical model (where averaging over the random fluctuations of the compound- 
nucleus resonances is an essential element), the question was addressed using 
a statistical model for the high-lying states. To this end, Feshbach's projec- 
tion operator method has been extended to the bound-state problem and 
used to estimate the resulting uncertainties in ground-state energies and wave 
functions [2]. (ii) A series of ever-refined nuclear mass formulas with a size- 
able number of parameters has been used to fit the known binding energies. 
In spite of many years of dedicated effort, there remains a difference between 
mass formulas and data: The data fluctuate about the (smooth) best fits. The 
fluctuations have a width of about 0.5 MeV and have been interpreted as being 
due to chaotic nuclear motion j3j. Since the spectral fluctuations of nuclear 
levels near neutron threshold (and near the Coulomb barrier for protons) follow 
random-matrix predictions T , the question arises whether this fact can be used 
to account for the observed fluctuations. 

We present a novel statistical approach to the question formulated above. 
In contrast to Refs. [2], it is not our aim to estimate the theoretical uncer- 
tainty of nuclear binding energies resulting from the use of a mean-field ap- 
proach. Rather, we are interested in the statistical fluctuations of ground-state 
and binding energies which are induced by those of high-lying conflgurations. 
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Figure 1: Schematic representation of the nuclear spectrum in the shell model. 
The semicircles represent the shells with n = I and n = 2. 



These represent genuine theoretical uncertainties which are quite independent 
of the particular approach used to calculate nuclear properties. To this end we 
describe the high-lying configurations in terms of a random-matrix model. In 
trying to be as conservative as possible, we only use the fact just mentioned 
that in many nuclei with mass number A > 30 or so, random-matrix theory 
correctly describes the statistical fluctuations of states at excitation energies of 
8 to 10 MeV. We show that from this input it is possible to estimate the mini- 
mum statistical uncertainties in nuclear ground-state energies and, from there, 
nuclear binding energies. We compare our results with the work of Refs. [2], 
and with the empirical results of Ref. jSj . 

Model. Starting point is the nuclear shell model. Neighboring major shells 
are separated by an energy difference Huj « MeV. We first neglect the 

mixing between major shells due to the residual interaction. Then, the nuclear 
ground state is obtained by diagonalizing the shell-model Hamiltonian (which 
consists of the single-particle energies of the lowest shell and of the matrix ele- 
ments of the residual interaction with respect to many-body states in that shell). 
The resulting ground state has energy Eq. The corresponding eigenfunction is 
composed of configurations belonging to the lowest major shell. We assume 
that the residual interaction within the states of each of the remaining major 
shells (excitation energies nhu) with n — 1,2,...) is so strong as to completely 
mix the configurations within that shell. Then, the spectra and eigenfunctions 
within each such shell are described by random-matrix theory, more precisely: 
by the Gaussian orthogonal ensemble (GOE) of random matrices. For n > 1, 
this assumption is consistent with the spectral fluctuations analyzed in Ref. 0], 
and with the properties of nuclear cross-section fluctuations at higher excita- 
tion energies The resulting nuclear spectrum is schematically depicted in 
Fig.Q] Each of the higher shells (n > 1) gives rise to a Gaussian-shaped mean 
level density. In what follows we replace each Gaussian by the semicircle char- 
acteristic of random-matrix theory. For our order-of-magnitude estimates, the 
difference should be irrelevant. We focus attention on the shell with n = 1 and 
later show that higher shells do not affect the result signiflcantly. We estimate 
the shift and the uncertainty of the ground-state energy Eg using (i) perturba- 
tion theory, (ii) the supersymmetry approach. Here, _Eo < is the difference 
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between the location of the ground state and the center of the shell with n = 1. 
The semicircle has radius 2A. Estimates for Eq and of A are given below. For 
simplicity, we focus attention entirely on the ground state and omit the other 
states in the shell with n — 0. 

Perturbation Theory. We denote by Voj the matrix elements of the residual 
interaction connecting the ground state with the states j — 1, ... in the 
major shell with n = 1, and by Ej the GOE eigenvalues in that shell, counted 
from the center. For N, the number of states in that shell, we have N ^ 1. 
The energy Eq of the perturbed ground state is given by 

Within the GOE, the Voj's and the Ej^s are uncorrelated random variables. 
The Vbj's are Gaussian distributed with zero mean value. Using the standard 
definition F-'^ = 27rV^^/d of the spreading width for mixing of the ground-state 
wave function with the states of the n = 1 shell, we write the second moment 
VojVok — (l/(27r))F^(i 6jk- Here, d = ttX/N is the mean level spacing in the 
n = 1 shell. Averaging Eq. over the GOE, we find to lowest order in T^, 

Eo^E, + ^J<lE'^-^j:^iE'-E,). (2) 



By definition, J^j^i^' ^ Ej) is equal to p{E'), the average level density nor- 
malized to N . Hence we find approximately 

fJ-a 

^^0 « i?o + ^ . (3) 

The level shift F''^A/(2i?o) is negative, as expected. To calculate the variance, 
we define p{Ei,E2) — J2jk ^(Ei — Ej)S{E2 — Ek) and use the fact that the 
correlation function p{Ei, E2) — p{Ei) p(i?2) tends to zero over distances \Ei — 
£^2! that are of order d (and not Nd). Proceeding as before, we find 
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[E,-E,Y-^E,-Eq^ =-^—^ . (4) 

Here c denotes a constant of order unity (and not of order N). We observe 
that with ^ 1, the fluctuations are much smaller than the level shift. Rough 
estimates are obtained by putting 2A w {2/i)fiU} (the semicircles of the shells n 
and n + 1 overlap), Eq « — (5/3)?iw, and F-l^ « 3 MeV. The last figure roughly 
corresponds to the spreading width of the giant dipole resonance which accounts 
for the mixing of one-particle one-hole and two-particle two-hole configurations 
and is, thus, similar to the mixing of the ground-state wave function with the 
wave functions in the shell with n = \ considered here. For the shift, this yields 
a typical value of —0.3 MeV. The r.m.s. fluctuation is roughly given by that 
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value divided by viV. To estimate N we use N = irX/dt- Here dt is the mean 
level spacing at neutron threshold. For heavy nuclei, dt « 10 eV, hco « 6.5 
MeV, and N « 10^ while for nuclei with mass number A around 50 we have 
dt « 10 keV, huj w 10 MeV and TV « 10^. We find that the r.m.s. fluctuation is 
about 0.3 keV for heavy nuclei and about 10 keV for nuclei with A w 50. These 
figures are order-of-magnitude estimates only, of course. 

Supersymmetry. We use the supersymmetry method (SUSY) [7| in order 
to show that the perturbative estimates given above essentially remain valid 
beyond the domain of perturbation theory. The calculations are rather technical 
and cannot be given here. A more detailed account is in preparation Suffice 
it to say that using SUSY and N ^ 1, we find the following saddle-point 
equation for the a matrix of the one-point function: 

^ A J_ A 

E-\<7^ N E-Xa~ {l/2)TiX/{E - Eq) ' ^ ' 

Here the first term on the r.h.s. is the usual contribution which, taken by 
itself, results in the semicircle law. (We recall that the spectrum is located in 
that energy interval where a has non-real values, with the average level density 
proportional to Im a). The second term is the correction due to the coupling 
with the ground state. This term yields an additional contribution to Im a in 
an energy interval of width AE located at Eq outside the semicircle and in the 
vicinity of Eq. We solve Eq. (O for the diagonal elements of the a matrix using 
perturbation theory. For the shift, we retrieve our previous result, Eq. 0. For 
AE we find 

. ^ 1 r^A 4A 

This agrees with Eq. Q except for the last factor. This factor is of order unity. 
It arises because in Eqs. Q and ©, we have calculated different expressions 
(the r.m.s. fluctuation and the width of the spectrum, respectively). We have 
used Eq. ISJ also to calculate shift and width non-perturbatively for values of 
Eq close to the edge of the semicircle. We have found that these values differ 
from the perturbative ones only by numerical factors. We conclude that our 
estimates in Eqs. Q and jSJ are approximately valid in a wide range of values 
for Eq and F^-. 

Discussion. So far, only the influence of the shell with n = 1 on the ground- 
state energy has been taken into account. What about higher shells {n > 
1)? Since SUSY has corroborated the perturbative results, we may safely use 
perturbation theory to answer the question. Formally, Eqs. Q and Q remain 
valid, with Eq replaced by the distance D„ between the energy of the ground 
state and the center of the shell n under consideration, and F-^ replaced by T^^ , 
the spreading width describing the mixing of the ground-state wave function 
with the configurations in shell n. While {Dn)~^ decreases only like n~^, F^ 
falls off very rapidly (exponentially). This is because F^ depends partly upon 
the overlap between the ground-state wave function and the configurations in 
shell n. In a different context, this exponential fall-off has been demonstrated 



4 



numerically in Figure 1 of Ref. jH] where the dependence of the ground-state 
energy upon the truncation of the shell-model space was studied. Therefore, 
we expect that the inclusion of higher shells would affect our estimates only by 
a factor two or three or so. It is possible, of course, that our treatment is too 
conservative, and that the domain where the GOE applies, extends further down 
in the spectrum. This would reduce \Eo\ and N and, thus, enhance both the 
shift and the r.m.s. fluctuation. Unfortunately, there is hardly any empirical 
evidence one way or another concerning that problem. In summary, we find 
that the shift of the ground-state energy is of the order of or larger than 0.5 to 
1 Mev for all nuclei, while the fluctuations amount at least to several 10 keV 
(1 keV) for medium-weight (heavy) nuclei, respectively. It goes without saying 
that these fluctuations do not describe actual fluctuations of the ground-state 
energy, but rather the minimum theoretical uncertainty due to lack of precise 
information about the states and interactions in the shells with n > 1. 

So far we have considered only a single level rather than the complete set 
of states belonging to the lowest shell. We do not expect that the presence of 
these states strongly affects our estimates. To be sure that this is indeed the 
case, we should also derive and solve the saddle-point equation for this case. 
This work is now under way [5]. 

How are our results related to those of Refs. The problem addressed in 
Refs. 121 was not to compute the r.m.s. fluctuations, but rather to relate the 
theoretical uncertainty of the binding energy per nucleon to that of the mean- 
field energy and to that of the residual interaction (the difference between the 
true interaction and the mean fleld) . The latter connects the mean-fleld ground 
state to the excited states. This was done in the framework of a somewhat 
simpler problem, the energy per particle of nuclear matter, rather than for 
finite nuclei. 

In most nuclei, the mean-field energy differs from the energy denoted by 
Eq in the present paper. Indeed, our definition of Eq includes the effects of 
the residual interaction within the lowest shell and differs, therefore, from the 
mean-field energy. Moreover, we focus attention on stochastic uncertainties 
resulting from the chaotic behavior of excitet states in the vicinity of neutron 
threshold. We do not attempt to estimate the error that comes with a theory 
based upon a mean-field approach. Therefore, the approach of Refs. 1' and the 
one used here, are rather different, and a direct comparison of the results is not 
possible. 

Fluctuations of Binding Energies. How can we use our results to estimate 
the ffuctuations in the binding energies of nuclei? The binding energy is the 
minimum excitation energy at which the nucleus fragments into A separate 
nucleons. We do not believe that it is possible at present to determine this 
excitation energy reliably in the framework of a nuclear-structure calculation 
for nucleus A. Indeed, to the best of our knowledge, nuclear-structure calcula- 
tions have rarely been used to calculate binding energies, except for the lightest 
nuclei. To answer the questions raised in Ref. P], we use the following indirect 
approach. Our calculations yield the fluctuation of the nuclear ground-state 
energy with respect to the center of the shell with n = 1. For most nuclei, the 
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threshold for neutron or for proton emission hes somewhere in that shell. Thus, 
we identify the calculated fluctuation with the fluctuation of the ground-state 
energy with respect to neutron (or proton) threshold. The fluctuation of the 
nuclear binding energy is then obtained by the combined effect of removing 
successively A nucleons from the nucleus. 

We apply this scheme in reverse order, starting from a nucleus with mass 
number and adding A — Aq nucleons to it. The first nucleus (in increasing 
order of A) where we have definitive evidence (both experimentally and theo- 
retically) that GOE statistics applies is ^^Al We assume that nuclei 
with mass numbers < 25 do not have any fluctuations in either ground-state 
or binding energy. Therefore, we choose Aq = 26 and an initial value 30 keV for 
the fluctuation of the ground-state energy and, thus, of the binding energy. To 
calculate the fluctuations of the binding energy for nuclei with mass numbers 
A > Aq we assume, for simplicity, that the fluctuations of the ground-state 
energy have a Gaussian distribution. (In the framework of our model, this may 
or may not be the case. We have not explored the solutions of the saddle-point 
equation to the point of ascertaining the answer to this question but recall that 
for the main part of the spectrum, the distribution has the shape of a semicircle 
and not of a Gaussian. Such details should, however, be irrelevant for a rough 
estimate). Under this assumption, the fluctuation of the binding energy of a 
nucleus with A nucleons is obtained by the convolution of (^ — ^o) Gaussian 
distributions, each with its own >l-dependent width. The resulting distribution 
is again Gaussian. Its variance (cta)^ is the sum of the variances of the indi- 
vidual Gaussians, the latter being given by the ^-dependent r.h.s. of Eq. Q). 
Thus, 

(..)'^,30,<cVf+ t (T) 

A'=Ao + l ^ 

Eq. (|7|) holds provided that the fluctuations in different nuclei are statistically 
uncorrelated. We are not aware of any evidence to the contrary. The most sig- 
nificant dependence upon A' on the r.h.s. of Eq. (0 resides in = nX/d. Since 
N increases strongly (nearly exponentially) with A, the individual terms under 
the sum decrease rapidly, and aA is a monotonically increasing function of A 
with ever decreasing slope. This statement applies in general, quite irrespective 
of any details of the A-dependence of N. We roughly model the A-dependence 
of the terms under the sum by writing 
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(30 keV)2 exp( - VA) ■ (8) 



This formula takes account of the initial condition at Aq = 26 and models the 
^-dependence in terms of the exponential which roughly reflects the depen- 
dence of the average level density on A. (A more precise expression would have 
exp ^/A replaced by exp[Q;\/^], with a of order unity and weakly dependent on 
A. We disregard such corrections in view of the rather rough overall approach 
we have taken). We use this expression in Eq. replace the summation by 
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a(A) 



O.f 



0.6^ 



0.4 ^ 
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Figure 2: Fluctuations of nuclear ground-state energies. Dashed curve: fit to 
the experimental data taken from Ref. '121; dot-dashed curve: formula (9) of 
the text; continuous curve: formula (9) with the factor 30 replaced by 120 keV . 



an integration, and obtain 



a{A) = (30 keV) 



i + 2(v/a;+i)(i- 



VA+1 



-exp{y^-VA}) 



1/2 



(9) 



In Figure|21 we compare our Eq. © with the result reported in Ref. |3] for the 

fluctuations of the binding energies versus A. To emphasize the uncertainties 
in our estimates, we show two curves, one of which represents Eq. JHJ as it 
stands while for the other, we have arbitrarily replaced the factor 30 keV on 
the r.h.s. of Eq. by 120 keV. We see that for heavy nuclei, the data and 
our result have similar values while for A « 50, the difference is larger. In 
view of the rather large uncertainties of our estimates, we would be satisfied 
with an order-of-magnitude agreement and do not attach much significance to 
such statements. More important, in our view, is the fact that we predict a 
monotonically increasing function of A while both, the data and the chaotic 
model of Ref. [HI, yield a monotonic decrease with A. As for the data, we can 
think of additional effects that might affect the ^-dependence. For instance, it 
is conceivable that for light nuclei, the clean separation between surface terms 
and volume terms assumed in the mass formulas, does not apply. This might 
add a correction which strongly decreases with A and which has to be removed 
from the data before a comparison with our results can be made. Within the 
conservative phenomenological approach which we have taken, there seems no 
way of escaping the conclusion that a a increases monotonically with A. This is 
clear from Eq. Q which, in turn, rests upon our identification of fluctuations 
of the binding energy of a nucleus with mass number A with the fluctuations 
of the ground-state energies with respect to neutron or proton threshold of all 
nuclei with mass numbers A' < A. 
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